A spin-dependent local moment approach to the 
Anderson impurity model 



Choong H Kim^ and Jaejun Yu^'^ 

^ Department of Physics & Astronomy and Center for Strongly Correlated 

Materials Research, Seoul National University, Seoul 151-747, Korea 

^ Center for Theoretical Physics, Seoul National University, Seoul 151-747, 

Korea 

E-mail: jyuOsnu.ac.kr 

Abstract. We present an extension of the local moment approach to the 
Anderson impurity model with spin-dependent hybridization. By employing the 
two self-energy description, as originally proposed by Logan and co-workers, we 
applied the symmetry restoration condition for the case with spin-dependent 
hybridization. Self-consistent ground states were determined through the 
variational minimization of ground state energy. The results obtained with of 
our spin-dependent local moment approach applied to a quantum dot system 
coupled to ferromagnetic leads are in good agreement with those obtained from 
previous work using numerical renormalization group calculations. 



PACS numbers: 73. 21. La, 75.30.Hr, 71.27.ta 
1. Introduction 

The Anderson impurity model (AIM) \X, and its extensions have been a matter of 
central importance in the recent developments of condensed matter physics. The AIM 
serves not only as a prototype model for the Kondo effect [2] but also as a key ingredient 
in the dynamical mean-field theory (DMFT) [31 H] for strongly correlated electron 
systems. Recent advances in quantum dot (QD) experiments have demonstrated that 
a quantum dot connected to leads can act as a magnetic impurity in metal so that 
the Kondo-type behavior emerges at low temperatures [S] . Indeed the Kondo effect in 
quantum dot systems have been probed by many theoretical [51 [7] and experimental 
studies [51 [21 [in] • The fine tunability of control parameters in quantum dot systems 
has spurred the investigation of Kondo physics in various aspects including a quantum 
dot coupled to ferromagnetic leads. A flood of very recent works [TT[ [T2l [T3l fT4l [TSl [T6] 
has focused on this issue as motivated by its potential applications to spintronics. 

Although there are many theoretical approaches to the solution of the AIM, it 
is still not easy to deal with the AIM coupled to ferromagnetic leads, where spin- 
dependent charge fluctuations need to be taken into account. In other words, it is 
necessary to get a non-perturbative method which can deal with both charge and spin 
excitation channels. For instance, some of slave-boson mean-field calculations could 
not describe the finite splitting of the Kondo peak [HI [T2] in the spin-dependent case 
due to the absence of charge fluctuations. Although the numerical renormalization 
group (NRG) calculation is known to provide accurate results for the impurity 
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Figure 1. Quantum dot coupled to ferromagnetic leads. 



problem, the standard NRG technique seems not good enough for the spin-polarized 
system due to the absence of energy scale separation [T7| . 

As a way to find an effective approach to the quantum dot coupled to 
ferromagnetic leads, we have considered the local moment approach (LMA) originally 
suggested by Logan and co-workers [TH [T^]. Based on the intuitive notion of 
local moment fluctuations and the symmetry restoration condition, this approach is 
technically simple and transparent and yet can cover all energy scales and interaction 
strengths ^U\. Since it relies on the degeneracy of the two mean-fleld saddle point 
configurations, however, it is not obvious how to handle the case with spin polarization. 
The original LMA of the Anderson impurity model in an external field [STJ [55] did not 
treat each spin component of the full Green's function separately, but just calculated 
the sum. In this work, we extended the LMA formalism to include the spin-dependent 
hybridization. The basic idea is a generalization of the two-self-energy description, 
allowing the variation of the weights of the spin-up and spin-down configurations. 
The weights of each local moment configuration can be determined through the 
minimization of the total energy of impurity system. 

This paper is organized as follows. In section [2l we give a general description of 
model Hamiltonian for the quantum dot coupled to ferromagnetic leads. In section [3l 
we describe the formalism of our spin-dependent LMA (to be abbreviated as sLMA) 
including basic ideas on the generalization of the two-self-energy description, the 
symmetry restoration condition, and the determination of configuration weights. The 
calculated spectral function for the asymmetric Anderson model is given in section [4l 
together with a comparison with previous works. 



Let us consider a model for the quantum dot (QD) coupled to two ferromagnetic leads. 
Figure [T] shows a schematic view of the QD system coupled to ferromagnetic leads 
with different spin polarization directions. We assume that a single level QD with a 
charging energy U. It is well known that the QD is equivalent to a single-level impurity 
state [5J [7] in the AIM. In this system, two cases with parallel (P) and antiparallel 
(AP) magnetic configurations are possible for the ferromagnetically ordered leads. By 
making the canonical transformation il3j , both P and AP configurations could be 
mapped onto an effective model with a single lead. After the transformation, the 
AP configuration was shown to be equivalent to a usual QD coupled to a single non- 
magnetic lead. Thus, from now on, we will focus on the P configuration where the 
spin-dependence becomes explicit. With the Fermi energy kept at the energy origin 
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{ep — 0), the Hamiltonian for the P configuration can be represented by 

ly ka 

where c\^ {ck<j) is the creation (annihilation) operator for an electron with the wave 
vector k and spin a in the leads. And d\ {d^r) is the creation (annihilation) operator 
for electrons in the dot. The local Zeeman coupling energy h — ^gfisH is included 
for the sake of generality. Here it is noted that the effective dot-lead coupling Vka is 
connected to the original QD model in figured] 

Vka = ^{VLka+VRka)- (2) 

Considering its parametric dependence, the spin-dependent tunneling amplitude V^k 
can be described by introducing a polarization parameter p for the spin-dependent 
host density-of-states (DOS) pa{(-o): 

P„{lj) = p{l + ap) for - L» < w < D. (3) 

Neglecting the dependence on k of Va-k in the large band width limit of conduction 
electrons, the spin-dependent hybridization parameter A^. = |T4P(cl'+ — e^^)^^ 
can be approximated by a simplified form: 

A,{io) ~ -iAo(l + (Tp)sgn(w) (4) 

where Aq — ttV^ p. Some authors have pointed out that the shape of the conduction 
band may be important in this problem |151 116j . However, for the sake of simplicity, 
we only deal with the Z? — > oo limit without considering the detailed band structure. 
Therefore, the Green's function of the d electrons can take a simple form of 

g^^{uj) ^uj+ - e, -ah + iAo(l + crp)sgn(w). (5) 

in the non- interacting limit (U — 0). In what follows, we consider the case oi h — 
only. 



3. Formalism 



3.1. Local moment approach 

The LMA starts off based on the mean-field, i.e., unrestricted-Hartree-Fock (UHF), 
solutions of the AIM. For large U values, the UHF treatment of the single-orbital AIM 
gives the doubly degenerate local moment solutions Ij. Two solutions are denoted 
by a = j4 or i?, corresponding to the local moments fi — +\fJ-\ and — respectively. 
In order to remedy the broken symmetry nature of the UHF solutions, LMA employs 
the two-self-energy description (TSE): 

G^{cu)^^[GAa{i0)+GBa{uj)]. (6) 

In this framework, the LMA is designed to describe the local moment at the impurity 
site, fluctuating between the two configurations, A and B. In addition, to include the 
dynamics of the self-energy going beyond the static UHF self-energy, the LMA uses 
the HF Green's function as a bare propagator for the calculations of the higher-order 
diagrams (which will be discussed in detail in section 13. 4p . The higher-order terms 
for the self-energy are indeed enough for describing the Fermi liquid behavior at low 
energy scale, but to do this it is required to introduce the symmetry restoration (SR) 
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condition and the Friedel sum rule for recovery of the correct fow energy behavior. The 
SR condition plays a role in linking the two configurations A and B at low energies. 
Taking advantage of its efficiency and transparency, the LMA has been applied to the 
lattice model within the framework of dynamical mean- field theory [231 1211 US [55] . 



3.2. Mean-field approximation 



In the absence of spin polarization, the UHF approximation gives two degenerate 
broken symmetry solutions: fi = -\-\^J-\ and — However, in the presence of 
spin polarization, the two mean-field saddle points become no more degenerate. To 
distinguish two solutions, we labeled two configurations hy a ~ A or B, corresponding 
to the local moment /i = +1/1^1 and — l/isj, respectively. 

Taking account of the configurations A and B, we obtained the mean-field 
propagator Q^a for a ~ A,B and cr =1, |: 



where 



Lu - CiA + + iAo(l -I- (Tp)sgn(w) 
(w) CJ+ - e,B - crxB + iAo(l -t- (Tp)sg\i{u) 



(7) 



(8) 



Here we introduced the configuration-dependent on-site energy, which may be 
necessary for describing the a-dependent renormalization of the on-site energy under 
the asymmetric configuration. And the corresponding spectral densities D^^{oj) are 

^L(^) = --sgn(Lj)Img„,(c^). (9) 

TT 

Consequently, the mean-field charge and moment could be determined from the self- 
consistent solution of the following equations: 

dujD^^^{uj\ eia,Xa) 



— oo 




dio 



where a 



\^la\ = a / duj Dl^{uj;eia,Xa) - D'l^{uj;eia,Xa) 

J — oo 

and — for A and B, respectively. 



(10) 



3.3. Two-self-energy description 

For the spin-dependent case, where the up-spin and down-spin symmetry no longer 
exists, we generalized the LMA description of the Green's function by allowing the 
two saddle point configurations to have different weights. Extending the original two- 
self-energy description, the average Green's function could be expressed as 

G^{lo) = CAGAaiio) + (1 - CA)GBaH, (H) 

where the Green's function for each configurations a is provided as 

G-i(c.)=g-jH-S„.H (12) 

with the a-dependent gaa and the corresponding self-energy Sckt(^): 

9aa{^) = w+ - e^a + iAo(l -\- crp)sgn(w). (13) 
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Figure 2. Principal contribution to the LMA Ect(lj). The shaded bubble 
represents the polarization propagator as given by the RPA-like particle-hole 
ladder sum in the transverse spin channel. 



The configuration weight ca or (1 — c^) reflects the probability of A or B configuration. 
In the absence of spin-dependence, ca is trivially equal to i. The dependence on a of 
gaa-i^) arises from the relative chemical potential shift of the A and B configurations 
due to the a-dependent renormalization of the bare energy level. 



3.4- Self- energies 

The evaluation of self-energy in the spin-dependent LMA (sLMA) followed the same 
steps prescribed in the standard LMA procedure, the details of which are given in 
[18j . The difference in our procedure lies mainly on the spin-dependent terms due to 
the non-degenerate local moment configurations. 

Following the notations in [18j, the self-energies were divided into the static HF 
term and the rest: 

^ ^{^a ~ acr\P-a\) + ^aa{i^)- (14) 

The purely static contribution is given by the UHF calculation. And all of the 
dynamics is contained in the I]q,o-(w) = Sqo-II^/qct}], which is a functional of the 
underlying MF propagator Gaait^), as shown in figure [D 



J-oo 27ri 



,iu; + u;i)UZi^i)- (15) 

This diagram contains the dynamical spin-flip processes, and within the random phase 
approximation (RPA) the polarization diagrams are calculated, 

1 - u °n„„(w) 

with the bare particle-hole bubble 

/oo 7 
-^gc.,{cJl)Gca{i0l-i0). (17) 
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3.5. Symmetry Restoration 

The symmetry restoration (SR) is required to recover the Fermi Hquid behavior at low 
CO. To achieve the same goal in the spin-dependent formalism, we have to modify the 
original SR condition for the generalized two-self-energy description with the different 
weights CA and (1 — ca) for spin-up and down configurations, respectively. In terms 
of self-energies, we can rewrite (jlip as 

1 CA ^ 1- CA 



(18) 



We could consider the following relation between the single self-energy and the 

two-self-energies T,aa{^^)'- 

ti -t- Scr(c^) = CA!Aa(y>) + (1 - CA)!Ba{'^) 

w+ + iAo(l + (Tp)sgn(w) - CAfsai^^) - (1 - ca)/a<tM ' 

where /acr(w) = eia + 'Saai^) is defined with a = A, B. Since the two-sclf-energies 
vanish at the Fermi level, i.e., l'mY,aa-{uj = 0) = 0, by imposing the condition that the 
imaginary part of self-energy at the Fermi level must be equal to zero, we could obtain 
the following generalized SR condition: 

e,A + = 0) = + Sg, = 0) . (20) 

This SR condition contains the relative chemical potential shift of two configurations. 
One can readily note that this generalized version of SR reduces to the original SR 
condition when there is no relative chemical potential shift, i.e., eiA = SiB- In addition, 
when the SR condition of ([^0]) is satisfied, the Fricdcl sum rule [Ul HH] equation can 
be written as 



t+Unag+^aai^ = 0;eia,Xa) (H-crp) Ao tan ^(1 - 2n; 



imp; 



.(21) 



This is complementary to the previous physical arguments related to the Fermi liquid 
behavior at the Fermi level. Here our self-energy arguments assert that the SR must 
be taken into account, when we are dealing with the two-self-energy descriptions. 

3.6. Filling constraint 

To make the practical calculation feasible, we introduced a physically motivated 
approximation, that the total impurity charge nimp does not depend on h or p. 

?^imp(^ = 0;p = 0) = ni,„p(ft,,p) = nf^^ih.p) = nf^^^{h,p) (22) 

where 







duj 



D^^{lo) + D^i{uj) . (23) 



Within this assumption, we fixed the value of rtimp as calculated for the h = p = 
system. Thus, we were able to determine CiA and CiB for a given n-i^p. 

Although the best way to do this is the self-consistent determination of n^^^p, there 
is no practical method available to determine the n^^^ self-consistently at present. 
Since the assumption on the independence of riimp on h and p is subtle but important 
in the real calculations, it may deserve more explanation. At least, both the exact 
Bethe ansatz solutions [29] and the NRG calculation [17] show the independence of 
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Figure 3. Energy vs. ca curve for {U/Aq = 12.5, ei/Ag = —5 and p = 0.0,0.1). 
The other configurations (e.g. p = 0.3, 0.5) has quite difli'erent energy scale to 
display in one figure. 



riiinp on h. Indeed this result supports that our assumption is vahd at least for the 
h ^ case. Even though the there is no explicit report for the case of p ^ 0, we 
can say that this assumptions is a good approximation for two cases where the charge 
fluctuations are strongly suppressed by either the particle-hole symmetry or the large 
on-site interaction, i.e., U 3> Aq. 

3. 7. Variational Principle 

For the ground state energy of the system, we adopted the following expression [50] : 



which is certainly valid only for the wide-band limit. It may be necessary to 
include the contribution from the conduction electron part for the general case. In 
the spin-dependent LMA, ca was taken as a variational parameter for the energy 
minimization. For a given ca, there exists a corresponding set of LMA parameters 
(ciA, CiA, Gis, xa, a;^) which satisfies the self-consistent equations (SR, Friedel 
sumrule and filling constraint). Hence we were able to calculate the ground state 
energy for the given ca which minimizes the total energy. 

Figure [3] shows the energy versus ca curve for U/Aq — 12.5 and ei/Ao = —5. 
(Here we display the results of ca > 0.5 considering the symmetry.) For the p — 
case, the minimum is present at ca ~ 0.5 as expected. It signifies that our energy 
variation scheme works for p — 0. For the case of p = 0.1, the energy minimum is at 
CA =i 0.77. 

4. Results 

We now turn to the sLMA calculation results. Figure H] shows the spectral functions 
7rAoD|(a;) and ■kAoDi{ui), i.e., local DOS of the QD, for the different values of 




(24) 
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Figure 4. Local DOS of the QD for the asymmetric Anderson model {U/Aq = 15 
and ei/Ao = —5). Note the different energy scales in two figures. 



lead polarization p in the P configuration. Please note the different energy scales 
in two figures. Our results on the finite splitting of Kondo peak are found to be 
consistent with those of numerical renormalization group (NRG) calculations [T3lll5j . 
As p increases, the Kondo peaks for ttAoD^{uj) and 7tAoDi{lu) shift in the opposite 
directions and the Kondo peaks split into two and the values of the spectral functions 
at Fermi level decrease. As a result, Kondo effects are reduced or suppressed in the 
presence of ferromagnetic leads. The relatively small magnitudes of the splitting are 
attributed to the similar underestimation of the width of Kondo peaks in the standard 
LMA. As shown in figure O qualitative features such as dependence on p of Kondo 
peak splitting S are in good agreement with NRG. One can observe that 6 is linear in 
the p, confirming the NRG result [13] . 

The finite splitting of Kondo peaks may be understood via Haldane's scaling 
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Figure 5. Splitting 5 of the Kondo peak as a function of p for U/Ag = 15 and 
Ei/Ao = —5. WK is defined as the FWHM of the p = Kondo peak. 



arguments: the renormalization of the on-site energy is spin-dependent due to the 
spin-dependence of the hybridization; 



where Ae was attributed to the sphtting of the renormahzed levels [15]. Then, the 
coupling acts as an effective magnetic field, leading to the finite splitting. 

In figure [H one can note that the mean-field peaks are also shifted in opposite 
directions. In the LMA scheme, the shift of the mean-field peak, arises from the 
different chemical potential shift of A and B, i.e., ^ eis- When we fixed the 
CiA = SiB = Ei, we could not observe the shifts although the weights of peaks become 
drastically different. The Hubbard satellites were more pronounced in sLMA, while the 
same mean-field peak shifts are rather small in the NRG calculations. It is well known 
that the high energy features are usually underrated in the calculation of dynamical 
properties with the NRG method ^T] . In addition, this problem becomes more serious 
in the spin symmetry breaking system due to the absence of energy scale separation 



In fact, the Kondo peak of the down-spin electron in our results is remarkable 
to observable. The peaks of down-spin is higher than the peaks of up-spin with the 
same p. Even though this is consistent with NRG results, there are some problems. 
In the NRG method, the height of the Kondo peak is decreased as increase of p, but 
our results show the opposite trend. The filling constraint which is introduced as the 
physically motivated approximation in section 13.61 could be the possible source of this 
difference. But it is very difficult to work out the source of the problem at this stage. 
Treating this problem could be an important issue for future work. 

5. Conclusion 

In this paper we introduced a spin-dependent local moment approach (sLMA) for 
the study of the Anderson impurity model with spin-dependent hybridization. As 
an extension of the standard local moment approach to the spin-dependent system, 
we employed the generalized two-self-energy description and symmetry restoration to 
deal with spin polarizations. The approach has been applied to a quantum dot system 
which is coupled to ferromagnetic leads. Our results for the asymmetric Anderson 
model are in a qualitative agreement with the NRG results. 




(25) 
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